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Abstract. Steady state fluctuation relations for dynamical systems are commonly derived 
under the assumption of some form of time-reversibility and of chaos. There are, however, 
cases in which they are observed to hold even if the usual notion of time reversal invariance 
is violated, e.g. for local fluctuations of Navier-Stokes systems. Here we construct and study 
analytically a simple non-smooth map in which the standard steady state fluctuation relation 
is valid, although the model violates the Anosov property of chaotic dynamical systems. 
Particularly, the time reversal operation is performed by a discontinuous involution, and the 
invariant measure is also discontinuous along the unstable manifolds. This further indicates 
that the validity of fluctuation relations for dynamical systems does not rely on particularly 
elaborate conditions, usually violated by systems of interest in physics. Indeed, even an 
irreversible map is proved to verify the steady state fluctuation relation. 
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1. Introduction 

One of the central aims of nonequilibrium statistical physics is to find a unifying principle in 
the description of nonequilibrium phenomena. Nonequilibrium fluctuations are expected to 
play a major role in this endeavor, since they are ubiquitous, are observable in small as well as 
in large systems, and a theory about them is gradually unfolding, cf. Refs. [|24l[Tn[T6l l9l[Tl[3l 
for recent reviews. A number of works have been devoted to the derivation and test of 
fluctuation relations (FRs), of different nature [[I8l[l9l|28l[3llll2ll2lil[34l. It is commonly 
believed that, although nonequilibrium phenomena concern a broad spectrum of seemingly 
unrelated problems, such as hydrodynamics and turbulence, biology, atmospheric physics, 
granular matter, nanotechnology, gravitational waves detection, etc. [[D |5l [6l |3, the theory 
underlying FRs rests on deeper grounds, common to the different fields of application. This 
view is supported by the finding that deterministic dynamics and stochastic processes of 
appropriate form obey apparently analogous FRs [[I] |3] 021 12 , and by the fact that tests of these 
FRs on systems which do not satisfy all the requirements of the corresponding proofs typically 
confirm their validity. Various works have been devoted to identify the minimal mathematical 
ingredients as well as the physical mechanisms underlying the validity of FRs Il29l l4ll30l[3ll. 
This way, the different nature of some of these, apparently identical but different, FRs has 
been clarified to a good extent flU [TH |2l |3l. However, analytically tractable examples are 
needed to clearly delimit the range of validity of FRs, and to further clarify their meaning. 

In this paper, the assumptions of time reversal invariance and of smoothness properties, 
required by certain derivations of FRs for deterministic dynamical systems, are investigated 
by means of simple models that are amenable to detailed mathematical analysis. In particular, 
we consider the steady state FR for the observable known as the phase space contraction 
rate A, which we call the A-FR, for dissipative and reversible dynamical systems, in cases in 
which A equals the so-called dissipation function f2 [[24|. and the A-FR then equals the steady 
state f2-FR BU. As will be shown below, the phase variables A and ^2 coincide provided that 
the probability density entering the definition of ^2 is taken uniform, as in the case of the 
equilibrium density for the baker map [|25l . Both the A-FR and the f2-FR rest on dynamical 
assumptions: While the steady state f2-FR has been proven to hold under the quite mild 
condition of decay of correlations with respect to the initial (absolutely continuous, with 
respect to the Lebesgue measure) phase space distribution dU, the A-FR has been proven 
for a special class of smooth, hyperbolic (Anosov) dynamical systems lfT8l[T9l , whose natural 
measure is an SRB measure. Indeed, there are almost no systems of physical interest that 
strictly obey such conditions. However, in a similar fashion, there are almost no systems 
of physical interest satisfying the Ergodic Hypothesis, and yet this hypothesis is commonly 
adopted and leads to correct predictions. Analogously to the ergodic condition, one may 
thus interpret the Anosov assumption as a practical tool to infer the physical properties of 
nonequilibrium systems. Nevertheless, it is important to investigate which aspects of the 
derivation of the A-FR are not essential to its validity. Along these lines one notices that 
the A-FR seems to inherently rely on a rigid notion of time reversibility, which, however, is 
not always satisfied Wf\ . and on the smoothness of the natural measure along the unstable 
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directions, which is also problematic. On the other hand, the validity of the f2-FR has never 
been explicitly checked on an exactly solvable model. 

By considering a fundamental class of chaotic dynamical systems, known as baker maps, 
we want to assess the relevance of the Anosov assumption and of time reversibility for the 
validity of the A-FR, in cases in which it coincides with the fi-FR. Also, by assessing the 
validity of these FRs while violating standard assumptions, we probe and extend their range 
of validity. The maps we consider are appealing, since they are among the very few dynamical 
systems which can be analytically investigated in full detail. For this reason, baker maps have 
often been used as paradigmatic models of systems that enjoy nonequilibrium steady states 
(NESS) [1251 [nmi [3511361. 

Our main results are summarized as follows: The assumed sufficient conditions of the 
standard derivation of the A-FR, i.e. smooth time reversal operator and the Anosov property, 
are not necessary. Indeed, the A-FR is verified in maps whose invariant measure and time 
reversal involution are discontinuous along the unstable direction. This result is connected 
with the fact that A equals Vt, and that the fi-FR is known to be quite a generic property of 
reversible dynamics. The Anosov condition allows the natural measure to be approximated 
in terms of unstable periodic orbits, which constitutes a convenient tool in low dimensional 
dynamics and even in some high dimensional cases liT4l [TTl l32l |5l. This approximation 
may hold even if the Anosov condition is not strictly verified, because periodic orbits enjoy 
particular symmetries which other trajectories do not [15] . However, if the Anosov condition 
is violated, one must check case by case whether the unstable periodic orbit expansion may 
be trusted. We will also face this issue, showing that in some case the unstable periodic orbit 
expansion becomes problematic, hence a different approach must be developed. In particular, 
we will profit from a separation of the full phase space into two regions, within each of which 
the invariant measure as well as the time reversal operator are smooth. Nevertheless, the full 
system is ergodic: the two regions are not separately invariant, and any typical trajectory 
densely explores both making the discontinuities relevant, e.g. for the role of periodic orbits. 

2. Time-reversibility for maps 

In this Section we review the concept of time reversibility for time discrete deterministic 
evolutions. In order to remain close to the notion of (microscopic) time reversibility of interest 
to physics, one usually calls time reversal invariant the maps whose phase space dynamics 
obeys a given symmetry. In particular, one commonly calls reversible a dynamical system 
if there exists an involution in phase space, which anticommutes with the evolution operator 

mm. 

In practice, consider a mapping M : U ^ U of the phase space U C , d E N, which 
evolves points according to the deterministic rule 

ai„+i = M(x„) , (1) 

where n is the discrete time. The set of points {x^, X2, x^, ■■■}, obtained by repeated 
application of the map M, constitutes the discrete analogue of a phase space trajectory of a 
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continuous time dynamical system and, indeed, each x„ could be interpreted as a snapshot 
of the states visited by a continuously evolving system. If M admits an inverse, Af^\ 
which evolves the states backward in time, like rewinding a movie, with inverted dynamics 
= M~^{x_^j^i), M is called reversible if there exists a transformation G of the phase space 
that obeys the relation 

GMG = ,GG = I, (2) 

where / is the identity mapping. 

This is not the only possible notion of reversibility; there exists a variety of weaker as 
well stronger properties [12211411 . which may be thought of as abstract counterparts of the time 
reversibility of the dynamics of the microscopic constituents of matter. For every x E U, the 
symmetry property Eq.© obviously implies 

GMGM{x) = x. (3) 

If M is a diffeomorphism, as often assumed [fTT|. Eq.(|3]) can be differentiated to obtain 

DG{MGM{x))DM{GM{x))DG{M{x))DM{x) = J , (4) 

where DM{x) denotes the Jacobian matrix of M evaluated at the point x of the phase space, 
and similarly DG{x). Using the relations [DM]~\x) = DM-\M{x)) and [DG]~\x) = 
DG{G{x)) leads to 

DM{GM{x))DG{M{x))DM{x) = DG{MGMG{x)) 



DG{M{x))DM{x) = DM-\MGM{x))DG{MGMG{x)) 



DM{x) = DG{GM{x))DM-^{MGM{x))DG{MGMG{x)) , (5) 
which, together with (|2l), yields 

DM{x) = DG{GM{x))DKr\G{x))DG{x) . (6) 
Moreover, computing the determinant of the matrices in Eq.® we obtain 

JM{GM{x))JM{x) ^''\^}f^ = 1 (7) 
Jg{x) 

where Jm{x) = \ det DM{x)\ and Jg{x) = \ det DG{x)\ stand for the local Jacobian 
determinants computed at x. Because the involution G is unitary and Jg{x) = 1 for every x, 
by definition, Eq.© can be simplified to obtain 

Jm{x) = J]:l{GM{x)) (8) 

for all X in the phase space. This equation provides a key ingredient for the derivation of 
fluctuation relations in dynamical systems |[T9l l20l ITTI . as we will also see later on for our 
examples. 

Reversible dissipative systems have been discussed extensively in connection with so- 
called thermostatting algorithms, both for time continuous |[8l [131 IH and time discrete 
[[2n [3511371 dynamics. Special attention has been paid in these systems to the time average of 
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the phase space contraction rate A(x) = — In Jm{x), which is an indicator of the dissipation 
rate. The other indicator recently used in connection with FRs is the dissipation function 
which, in our context, takes the form [31 [l6l 

for a given phase space probability density p. Obviously, Cl takes different forms depending 
on p, and one has A = f2 if p is uniform in the phase space, which will be our case. Hence, in 
the following we only use A for simplicity. 

On a trajectory segment of duration n steps, starting at initial condition Xq, the time 
average of A is defined by 

^ 71—1 

K{x^)--J2^nJM{M''{xo)) ■ (10) 

fc=0 

Given this trajectory segment, let us call reversed trajectory segment the segment of duration n 
and initial condition GM'^(xq) = M~^G{xj^), cf. Eq.©. Its average phase space contraction 
rate may be written as 

^ 71—1 

A„(GM"(xo)) = --Y,^nJM{M^GM^{x^)) 

k=0 



-- VlnJM(GM-'=+"(xo)) 
n ^-^ 

k=0 

^ 71 — 1 

-5^1nJM(M-'=+"-i(xo)) (11) 



n 

k=0 



in which the last equality follows from Eq.® if the dynamics is time reversal invariant. We 
have thus shown that the phase space contraction rates of reverse trajectories take opposite 
values, 

A„(GM"(xo)) = -A„(xo) , (12) 

in time-reversible dissipative systems. It is interesting to note that, in discrete time, the initial 
condition of the reverse trajectory is constructed by applying the reversal operator G to a 
point, M"(xq), which is not part of the forward trajectory segment, but is reached one time 
step after the last point of the original segment. This equation is at the heart of the proof of 
steady state FRs for reversible dynamical systems. 



3. The A-FR 

The steady state A-FR was first obtained by Evans, Cohen and Morriss [|T8l for a Gaussian 
ergostatted (i.e. constant energy [|8l) particle system, whose entropy production rate is 
proportional to the phase space contraction rate. It was then rigorously shown to be 
characteristic of the phase space contraction rate of time reversal invariant, dissipative, 
transitive Anosov systems by Gallavotti and Cohen [[T9l . 
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This relation may be expressed as follows. Consider the dimensionless phase space 
contraction rate, averaged over a trajectory segment of duration n, with middle point x, in the 
phase space U, 

n/2-l 

where, without loss of generality, n is even and 



(A) = / A(x) /i(dx) 

J M 

is the nonequilibrium steady state phase space average of A, computed with respect to the 
natural measure /i on U, i.e. the M-invariant measure characterizing the time statistics of 
trajectories typical with respect to the Lebesgue measure. Then the Fluctuation Theorem may 
be stated as follows lITQlfTTl : 

Gallavotti-Cohen Fluctuation Theorem. Let M be a C^^"^, a > 0, reversible Anosov 
dijfeomorphism of the compact connected manifoldU, with an involution G and a G-invariant 
Riemann metric. Let fi be the corresponding SRB measure, and assume that (A) > with 
respect to [i. Then there exists p* > such that 

p-5 < hm — — log— y —jr^<p + 5 (14) 

if\p\ < p* and 6 > 0. 

Eq.(fT4l). usually considered for an arbitrarily small 6 and by specifically dealing with 
the phase space contraction rate as an observable, refers to what we denoted as the A-FR 
in the introduction. According to this terminology, one may say that the Gallavotti-Cohen 
Fluctuation Theorem proves the A-FR under specific conditions. This theorem is a rather 
sophisticated result, obtained by heavily relying on properties of Anosov diffeomorphisms, 
hence, in principle, it is hardly generic (see also [|20l ). For instance, Ruelle's derivation [[TT| 
makes use of Bowen's shadowing property, topologically mixing specifications, properties 
of sums for Holder continuous functions, expansiveness of the dynamics, continuity of the 
tangent bundle splitting, the unstable periodic orbit expansion of fi, and large deviations 
results for one dimensional systems with short range interactions. In these derivations, time 
reversibility and transitivity are necessary to ensure that the denominator of the fraction in the 
A-FR does not vanish when the numerator does not, while the smoothness of the invariant 
measure along the unstable directions, which allows the periodic orbit expansion, is included 
in the SRB property of /i. Recently, Porta has shown for perturbed cat maps that the A-FR 
requires the existence of a smooth involution representing the time reversal operator. 

Experimental and numerical verifications of relations looking like Eq. (fT4)) . for 
observables of interest in physics, have been obtained for systems which may hardly be 
considered Anosov [[24l[T6ll23l . Therefore, especially in view of the fact that the observable 
of interest is not A, except in very special situations, various studies have argued that strong 
dynamical properties, such as those required by the standard proof of the fluctuation theorem 
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for A, should not be strictly necessary [|26l |4l |T6l [Tl. Indeed, according to these references, 
time reversibility seems to be the fundamental ingredient for fluctuation relations of the 
physically interesting dissipation, since a minimum degree of chaos, such that correlations 
do not persist in time, can be taken for granted in most particle systemsJI 

Here we will proceed to show that properties implied by the Anosov condition, like the 
smoothness of the natural measure along the unstable directions, are violated in some simple 
models while the A-FR still holds. 

4. The A-FR for a simple dissipative baker map 

Research on chaos and transport has strongly benefitted from the study of simple dynamical 
systems such as baker maps [|25l [T2l l35l HJ. These paradigmatic models provide the big 
advantage that they can still be solved analytically, because they are piecewise linear, yet they 
exhibit non-trivial dynamics which is chaotic in the sense of displaying positive Lyapunov 
exponents. There are two fundamentally different ways to generate nonequilibrium steady 
states for such systems flU, namely by considering area preserving, 'Hamiltonian-like' maps 
under suitable nonequilibrium boundary conditions [fT2l |38l [371, or by including dissipation 
such that (A) > 0, as required by the A-FR [[35ll2Tl l25l . Within the framework of the former 
approach, FRs for baker maps have been derived in Refs. [|35l l36l . Here we follow the latter 
approach by endowing the map with a bias, which can be represented by a suitable asymmetry 
in the evolution equation. This bias may mimic an extemal field acting on the particles of a 
given physical system by generating a current One should further require the map to 
be area contracting (expanding) in the direction parallel (opposite) to the bias, which is the 
situation in standard thermostatted particle systems [[211 [35l . 

We now discuss the proof of the A-FR for maps of this type. The probably most simple 
model is described in Refs. [[39l l37l l25l . Here we give, in a different fashion than in the 
book Ref. Il25ll . the proof of the A-FR for this system by including the one sketched in this 
book. This sets the scene for a slightly more complicated model, which we will analyze 
in the following section. The calculations that follow allow us in particular to investigate 
the applicability of the unstable periodic orbit expansion for cases in which the smoothness 
conditions that guarantee their applicability are violated, but to a different extent in the 
different models. 

Let U = [0, 1] X [0, 1] be the phase space, and consider the evolution equation 




ryn J ' 
{xn - l)/r 
r + lyn 



) 



for < a; < /; 



for / < X < 1. 



(15) 



X Of course, one may expect exceptions to this rule in cases where randomness in the dynamics is somewhat 
suppressed 141 [T6l . 
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At each iteration, U is mapped onto itself, and the Jacobian determinant is given by 

J A = f/^i for < X < /; 



Jm{x) 



= l/r = J-\ for / < X < 1. 



(16) 



The map M is locally either phase space contracting or expanding. Furthermore, the 
constraint r + / = 1 makes the map reversible, in the sense of admitting the following 
involution G, meant to mimic the time reversal invariant nature of the equations of motion 
of a particle system. 



Vg 



G 



X 

y 



1 — X 



(17) 



The map G amounts to a simple mirror symmetry operation with respect to the diagonal 
represented in FigH] 



1 



A 



B 



1 



G 



\-l 



I 1 

Figure 1. Involution G for the map defined by Eq.(fT5Tl. 
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Figure 2. Check of reversibility for the map in Eq.(fT5t. performed by verifying Eq.(|2|l. 
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The relation Jb = J^^ in Eq.(fT6l) is a direct consequence of the time reversibility of the 
model. To see how this occurs, let us first observe, with the aid of Fig. [2l that the following 
relations hold for the map Eq.(fT5Tl: 

GMA = B , GMB = A . (18) 

Combining this with Eq.®, we immediately obtain Eq.(fT6l). Relation ([8]) can be further 
exploited by introducing the Jacobians of the dynamics restricted to the stable and unstable 
manifolds in the generic regions i = {A, B}, which we denote by and J^, respectively. 
One then has 

which, considering the specific constraints of ourmap, {Ja)~^~^{Jb)~^ = 1 ^i^d Ja + Jb = 1' 
leads to 

JX = {JlV and J% = {Jl)-' . (19) 
These equations constitute a consequence, like Eq.®, of the time-reversibility of the model 

ED. 

A probability density p„ on W, given at time n, evolves according to the Frobenius-Perron 
equation as ll25l[T0ll 

Qn+l (M (X) ) = J^-/ (X) Qn (x) • (20) 

Correspondingly, the mean values of a phase function O : U ^ IR evolve and can be 
computed as 

{0)n = [ 0{x)dfXn{x) [ 0{x)pn{x) dx . (21) 

Ju Ju 

If {0)n converges exponentially to a given steady state value {O), for all phase variables c|§| 
one says that the state represented by the regular measure corresponding to the density p„ 
converges to a steady state, which yields the asymptotic time statistics of the dynamics. This 
state will be characterized by an invariant measure p, which typically is a natural one. For our 
models this measure is singular, because M is dissipative [|39l[25l . 

However, due to the definition of the map Eq.(fT5l). which stretches distances in the 
horizontal direction —the direction of the unstable manifolds—, every application of the map 
smoothes any initial probability density in that direction, so that our invariant measure is 
uniform along the x-axis. Therefore, to compute steady state averages it is not necessary to 
use the full information provided by the n — )• oo limit of Eq.(l20l). Without loss of generality, 
we may assume that the initial state is "microcanonical", i.e. its density is uniform in U, 
Po{x,y) = 1. Then each iteration of the map keeps the density uniform along x, while 
it produces discontinuities in the y direction, so that the n-th iterate of the density can be 
factorized as 

Pn{x,y) = C ■ pn{y) , (22) 

§ The space of phase functions depends on the purpose one has in mind. The choice of Holder continuous 
functions is common fTTl. 
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where p„ is a piecewise constant function, which gradually builds up to a fractal structure, and 
C is a constant that is easily computed to be 1 by requiring the normalization of p„. Hence, 
the varying averages of observables are computed as 

f dx f dyO{x,y)pn{y) , (23) 
Jo Jo 

and their steady state values are obtained by taking the limit n — )■ oo. The average of the 
phase space contraction rate, which is constant along the y-axis, is then easily obtained as 

(A)„ = - dx dy pniy)lnJ{x) 



= I dxln^+ I dxln-(/ - r)ln(//r) . (24) 
Jo I Ji r 

As this result does not depend on n, it does not change by taking the limit, and we have 

(A) = (/ — r) ln(//r), which vanishes for / = 1/2 and is positive for all other / G (0,1). From 

Eqs. (fTOl) and (fT6l) . we can write 

nAn = {a-f3) In Jb , (25) 

where a and (3 = n — a denote the number of times the trajectory falls in region A or region 
B, respectively. 

To proceed with the derivation of the A-FR for this map, one may now follow two 
equivalent approaches. First of all, observe that our map is of Anosov type, except for an 
inessential line of discontinuity, which does not prevent the existence of a Markov partition. 
Therefore, two basic approaches to the proof of the A-FR may be considered: One may either 
trust the expansion of the invariant measure in terms of unstable periodic orbits [|T4l ITSl . or 
one may adopt a stochastic approach to the fluctuation relation f25], motivated by the fact 
that our baker map is isomorphic to a Bernoulli shift, i.e. to a Markov chain whose transition 
probabilities fulfill 

where ^^'"(xo)' with A; G [0, n — 1], denotes the region containing the point M'^(xq), out of the 
two regions {A, B},p{ij^[k(^,j.^y, k — )■ i^k+if^^^y, k+1) denotes the probability that the evolution 
touches region «A/fe+i(2,.^) at the time step k + 1, given that it visited the region iMf^{x^^) the 
previous time step k, and p{iM''+^(xj^)'i k + 1) the probability that M'^+^(xo) belongs to the 
region ^m'^+H^^o)" the n — > oo limit, the latter becomes the invariant measure Ati^^fe+i^ ^ of 
the region iMf'+^ixj^) itself. 

If one uses unstable periodic orbits, the argument proceeds as follows: Every orbit uj is 
assigned a weight proportional to the inverse of the Jacobian determinant of the dynamics 
restricted to its unstable manifold, which is = {Ja)"{JbY ^ if ^ f^U^ region A a 
number a of times and falls in region B a number (3 of times. Then the probability that 
the dimensionless phase space contraction rate e„, computed over a segment of a typical 
trajectory, falls in the interval Bp s = {p — 6,p + 5), coincides, in the large n limit, with the 
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sum of the weights of the periodic orbits whose mean phase space contraction rate falls in 
,5. Denoting this steady state probability by 7r„(i?p 5), one can write 

vr„(i?p,5)^^ W > (26) 

where Nn is a normalization constant, and the approximate equality becomes exact when 
n -7- 00. Because the support of the invariant measure is the whole phase space U, time 
reversibility guarantees that the support of 7r„ is symmetric around 0, and one can consider 
the ratio 

T^njBp^s) ^ Yju,e,,{w)(iBp^s{J^) ^^^^ 

where each to in the numerator has a counterpart in the denominator, and the two are related 
through the involution G, as implied by Eq.(fT2l). Therefore, considering each pair of trajectory 
segments cu and ZJ, of initial conditions Xq and G'M"(xq) respectively, Eqs.(fT2l) and (fT9l ) imply 

e„(a;) = -enip) , (J^)-' = , (28) 

where, for sake of simplicity, by e„(a;) we mean the average of e„(xQ), based on any point 
Xq of the orbit u. Consequently, exponentiating the definition of e„A„/ (A), and recalling that 
JujJ^J^^ for every orbit cu, we may write 

f 7^ = f «-P[-((A)j> + a (29) 

where |e^| < 5 if e„(a;) G Bp g. Because each forward orbit uj in the denominator of Eq.(l27]) 
has a counterpart to in the denominator, and Eq.(l29l) holds for each such pair, apart from an 
error bounded by S, the whole expression Eq.dTTll takes the same value as each of the ratios 
Eq.(l29l). with an error |e| < 5, 

7r„(-Bp,5) _ ^„((A)p+e) 



e-^^^^>^^^>, (30) 

T^nyB-p^S) 

where e can be made arbitrarily small by taking 5 sufficiently small and n sufficiently large. 
For a given S, n must also be large because, at every finite n, the values which Cn takes 
constitute 2n+ 1 isolated points in [— 1 , 1] . Therefore, 7r.„ {Bp^s) vanishes if none of these values 
falls in Bp s, making the expression senseless. But the set of these values becomes denser and 
denser as n increases. Taking the logarithm of Eq.d30]). for consistency with Eq. (fT4l) . and 
choosing p among the values e„ which may be attained along a periodic orbit of period n, we 
may now write 

' ln^^#4 = en4v(«-/3)lnf-l OD 



n(A) 7r„(5_e„,5) "n(A)' \r ^ 

for any 5 > 0. The n — )■ 00 limit of the above expressions confirms the validity of the A-FR, 
under the assumption that the unstable periodic orbit expansion could be applied. 

From the point of view of the Bernoulli shift description we obtain the same result, 
supporting the applicability of the the unstable periodic orbit expansion, despite the 
discontinuity of the dynamics at /. Indeed, observe that / equals the probability = j\ yu(dai) 
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Figure 3. Illustration of the generalized baker map. Green lines: Piecewise linear one- 
dimensional map, which generates the dynamics along the unstable manifold. 



that the trajectory can be found in region A, and r equals the probability jiB 
it is found in region B. Therefore, one may write as well 



In 



In 



a 13 



P a 



Ib /^('i^i) that 



(32) 



which is due to the instantaneous decay of correlations in the Bernoulli process. This leads us 
to conclude that the violation of the Anosov property, in this simple baker model, is irrelevant 
for its behavior. 



5. The A-FR for a generalized dissipative baker map 



We now propose a novel, generalized baker map, which is different from previous models 
||25l |2T1 1371 [391 |35l by generating a discontinuity in the invariant density along the x-axis. 
As illustrated in Fig. |3l this is achieved by the map acting differently on four subregions of 
U = [0,1] X [0,1], defined by 



for < X < /; 






for / < X < |; 



for I < X < |; 



for I < X < 1. 



.(33) 



In the sequel, unless stated otherwise, by M we refer to the map introduced in Eq.(|33l). 
The model is fixed by choosing the value of / G [0, |], i.e. the width of sub-region A. 
The parameter which determines the dissipation, hence the nonequilibrium steady state. 
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corresponds to a bias h which is suitably defined by 6 = — J]^ = 2 — According to 



its geometric construction shown in Fig. |3] the map Eq.(l33T) is area-contracting in region B, 
area-expanding in region C, and area preserving in regions A and D. This is confirmed by 
computing the local Jacobian determinants of the map to 



Ja 
Jb 
Jc 
Jd 



[2(1-2/)]- 

2(1-2/), 

1, 



for 
for 
for 
for 



< X < /; 

1 <x <\; 

\<X<1. 



(34) 



The following involution G, 





2 

2x 
y 

2 

2x 



for 



< x < i; 



(35) 



for j: < X < 1. 



constitutes a time reversal operator for the map M defined on the unit cell. It consists of the 
composition G = F o S of two other involutions, with 5* permuting the left and the right halfs 
of the unit square, and F mirroring the regions along their respective diagonals for all values 
h G (-oo,l],cf.Fig.|4l 

G^FoS 



A B 




► CD 



A B 



1/2 



/ 1/2 3/4 1 1/4 1/21/2 + / 1 

Figure 4. Involution G for the map defined by Eq. 



c 


A 


B 


D 



1-2J 



1/2 



Analogously to Eq.(fT8l) for the map Eq.(fT5l). for the generalized map Eq.(l33l). Eq. (|35l) 
entails the relations 

GMA = A , GMD = D , GMB = G , GMG = B , (36) 

which can also be inferred graphically from Fig. |5] It is readily seen, again, that the Jacobian 
rule Eq.®, supplemented by Eq.(l36]). implies the relations Eq.(l34l). 

Let us now lift this biased dissipative baker map onto the whole real line in form of 
a so-called multibaker map, which consists of an infinitely long chain of baker unit cells 
deterministically coupled with each other. Multibakers have been studied extensively over the 
past two decades as simple models of chaotic transport [[T2l |9l |35l [2T1 [271 [37]|. In our model, 
which we denote by Mmb, all unit cells are coupled by shifting the regions B and G to the, 
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Figure 5. Check of reversibility for the map of Eq.(l3Jt. performed by verifying Eq.©. 















c 


D 











mh 




Figure 6. Illustration of a multibaker chain based on the unit cell defined in Fig. [3] featuring a 
flow of particles from the regions B and C of the cell m into, respectively, the neighboring cell 
TO + 1 on the right and onto the neighboring cell m — 1 on the left. The net flow of particles 
corresponds to the current ^P, which is found to be proportional to the average phase space 
contraction rate (A). 



respectively, right and left neighboring cells, cf. Fig. |6l Choosing 6 7^ 0, i.e. / 7^ \ then 
implies the existence of a current defined by the net flow of points from cell to cell. The 
map Mmb is area contracting (expanding) in the direction (opposite to the direction) of the 
current, analogously to the case of typical thermostatted particle systems [1211 This can 
be inferred from the graphical construction in Fig. ^complemented by the relations Eq.(l34l). 

To asses the validity of the A-FR for this model, let us observe that the form of the 
invariant probability distribution along the y-direction (the direction of the stable manifolds) 
is irrelevant, analogously to the case discussed in SecJH because the phase space contraction 
per time step. A, does not depend on y. By introducing the shorthand notation = In Jc we 

II Differently, the pump model of Ref. ll33l may be tuned to expand phase space volumes in the direction of the 
current. 
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have 



A(x) = A(a;, y) < 



for < X < /; 
for I < X < 



1. 



2' (37) 

-0, fori<a:<|; ■ 

0, for I < X < 1. 

The y-coordinate may then be integrated out, and one only needs to consider the projection of 
the invariant measure on the x axis, the direction of the unstable manifolds, which has density 

Px- 

The calculation of this invariant density can be conveniently performed by introducing 
a Markov partition of the unit interval, which separates the region < x < 1/2 from the 
region 1/2 < x < 1. Denote by pi and pr the projected density computed in these two 
regions and let T be the transfer operator associated with the Markov partition. One may then 
compute the evolution of the projected densities, which are now piecewise constant, if the 
initial distribution is uniform on the unit square. In this case the corresponding Frobenius- 
Perron equation Eq.(l20t takes the form [[T2l [ 

Pi{x^n+i) \^f Pi{x T= \ ^ r T 1 • (38) 

Pr{Xn+l) J \ Pr[Xn) 

According to the Perron-Frobenius theorem, the transfer matrix T has largest eigenvalue 
A = 1, whose corresponding eigenvector yields the invariant density of the system to 




p{x) i i+f - ^ - - 2' _ (39) 

^ Prix) = tI^, for 1 < x < 1. 

This result confirms that, by construction of the model, and differently from the one 
considered in Section |4l the density of the map Eq.(|33]) is not uniform along the x-direction, 
that is, it is actually discontinuous along the unstable direction. 

By using this density, the average phase space contraction rate can be calculated to 

(A) = -$(6)ln^>0 , (40) 

where 

^(6) = — ^ (41) 
is the steady state current in the corresponding multibaker chain. Note that 

*(&)^^ (b^O) , (42) 
hence we have linear response and a caricature of Ohm's law. Accordingly, we get 

(A)^^ (6^0) (43) 

for the average phase space contraction rate, as one would expect from nonequilibrium 
thermodynamics if this quantity was identified with the nonequilibrium entropy production 
rate of a system [|2T1 l35l . This confirms that our abstract map represents a 'reasonably 
good toy model' in capturing some properties as they are expected to hold for ordinary 



for < X < h' 
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nonequilibrium processes. Related biased one-dimensional maps have been studied in 
Refs. [SOlllll. Note that ^ = for / = 1/4, respectively 6 = 0, only, in which case the 
dynamics is conservative, and the model boils down to a special case of the multibaker map 
analyzed in Ref. [l38l. 

In order to check the A-FR for this model, we first need to define the transition 
probabilities pij of jumping from region i to region j, with i,j E {A, B,C, D} denoting 
the finite state space. They constitute the elements of the transition matrix 



/ 1 i\ 



P 



(44) 



21 1-21 
1 i 

\2i 1-21 y 

Note that P defines a stochastic transition matrix, which acts onto vectors whose 
elements are the probabilities to be in the different regions, in contrast to the topological 
transition matrix Eg. (1381). which acts upon probability density vectors. The left eigenvector of 
P, associated with the eigenvalue 1, corresponds to the vector of the invariant probabilities /ij 
of the regions A, B, C and D. Alternatively, since the projected invariant probability density 
is constant in each of these four regions, the /i/s are also immediately obtained by multiplying 
the relevant invariant density Eq.(l39]| with the width of the respective region. One way or the 
other, we obtain 

f ^, if i = A,C,D: 

The discontinuity of the invariant density Eq.(l39l) along the unstable direction, for / ^ 1/4, 
means that the Anosov property is more substantially violated here than for the map in Section 
m Therefore, the periodic orbit expansion used in Section |4] cannot be immediately trusted, 
and an altemative method is better suited to prove the validity of the A-FR. 

We may begin by considering a trajectory segment of n steps, which starts at Xq E and 
ends in x„ E ix^, hence visits the regions {ix^, ...,ix^}. Consider the first (n — 1) transitions, 
corresponding to the symbol sequence {ix^, and treat separately the last transition 

^ Denote by riij the number of transitions from region i to region j, along the 
trajectory segment of (n — 1) steps, and by Ui = J^ij-pi^j^o} "^v ^''^^^^ number of transitions 
starting in i. Some transitions are forbidden, as shown by Eq.(l44]|. hence the following holds: 

Uac + nAD + riBA + nBB + ricc + ncp + npA + ^gi? = n-l . (46) 

riA riB nc nu 

It also proves convenient to introduce the following symbols: 

.0, if the trajectory does not start in i; ^^^^ 
^* 1, if the trajectory starts in 2. 

, 0, if the trajectory does not end in i; , 

^^-^^ 1 -f.u . . A ■ ■ (48) 
1, if the trajectory ends m z. 
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and Ajj = n^i — nj^. The quantities n^j and rij^ take into account the possibility that the 
trajectory segment may, respectively, start in, or end into, the region i. Thus, we may write 
the following /?MJc balances 

J2 = - Ajj, Vj (49) 

for each region of the map. Next, we introduce the quantity 

g = ub - nc + ub-^ - nc-^ , (50) 
which lies in the interval [—n, n] and is related by 

Tin = g<P/n (51) 

to the average phase space contraction in a trajectory segment of n steps. 

To evaluate the ratio of probabilities appearing in the A-FR, let us denote by the 
region containing the point x, out of the four regions {A, B, C, D}, and let us focus on a 
single trajectory of initial condition Xq = {xo,yo) E i^^. For a given n, the sequence of 
transitions which take this point from region i^^ to region iM{xg) = ix-^, from region to 
region in^^^) = i^^ and eventually from region ZAf"-i(xo) = ^3^n-i region 2m"(xo) = 
does not depend on yo. The larger the value of n, the narrower the width of the set of initial 
conditions whose trajectories undergo the same sequence of n transitions experienced by the 
trajectory starting in Xq. Let u{xjj,n) = {x E U : M^{x) E iM'={xg)^f^ = 0, ...,n} C ix^ 
denote this set of initial conditions. The expansiveness of the map implies 

lim u{xQ,n) = {x = {x,y) : x = Xo,y E [0, 1]} . 

n— ^oo 

Because the phase space contraction A(x^) only depends on the region i^^ from which the 
transition i^^ — )■ i^^^^ occurs, all trajectory segments of n steps originating in uj{xq, n) 
enjoy the same average phase space contraction A„. The amount A„ is also produced by 
the trajectory segments which visit the regions ix^, %„_i ^nd eventually land in i^^, where 
ix^ 7^ is the other region reachable from ix^_^- Let Ci;(xq, n) be this second set of initial 
conditions producing A„ in n steps. The point Xq lies in i^^, i.e. i^^ = i^^, but differs from Xq 
and does not belong to uj{xq, n). Denoting by vr^(xg,n) the invariant measure of oj{xq, n), one 
finds: 

'^u){xQ,n) 

n-2 

= /^i^o JJp(^M%; k iM^+^x^; k + l)p{iM"-i{x^),n- 1 iM"(xo);n) 

fc=0 



,,. nBA+noA riBB+riDB riAC+nccnAD+ncD (■ , ■ „ — 1 , ^■ 

l^ts.oPDA Pbb Pcc Pad P[iM"-^{xf,),n - i ^ iM"{xf,), 



n) 



,, ^'"■A-^AA^nB-^BB^nc-^CC ^no-^DD w 

H'ixoPDA Pbb Pcc Pad ^ 



X p(iM"-i(xo); n-1 Hf"(xo); n) , (52) 
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where we made use of Eqs.(|49]l and of the equalities pij = pkj for all % ^ which can be 
deduced from an inspection of Eq.(l44]|. Similarly, one has 

7r^(xo,n) = l^i^/'DA^^^Vl^B^''''Vcc^''°V'AD^''°V{^M^-^^^^ -> -(53) 

Given the similarity of the expressions Eqs.(l52]) and (|53T ). and the fact that 

- 1 iM^{x^)\n) + PiiM'^-HxoY^n - 1 z^^; n) = 1 , (54) 
it is convenient to consider the set 

Ua;(xo,n) = cj(xo,n - 1) (55) 
whose measure is given by 

riA-^AA riB-ABB nc-Acc riD-AoD (Cf.\ 
7rw(Z(,,n-l) — H'i^PdA PbB PcC PaD ■ W^^^ 

This measure represents the contribution to the probability of producing A„ in n steps, given 
by the trajectory segments whose initial conditions lie in u:{xi^,n — 1). The steady state 
probability of A„ is then the sum of contributions like Eq.d56]). for all remaining sets of 
trajectories compatible with A„, characterized by distinct sequences of n — 1 transitions. 

As we discussed at the end of Section 2, for any initial point Xq in the phase 
space that experiences a mean phase space contraction A„(xo) in n steps, the point 
Xq^ = GMM"~^(xq) = GM^{x^) experiences the opposite mean phase space contraction 
A„(xq^) = — A„ (xo), cf. Eq.(fT2l). The trajectory segment of n steps, starting at Xq^j, is thus 
the time reversal of the one starting at Xq, and ^(^^^^(xq), is the set of initial conditions 
of the time reversals of the segments beginning in uj{xj^,n). The segments beginning in 
a;(GM"(xo),n) visit the regions iau^ix^), icM^-Hx,-,), ■■■ ^Gfeo). hence they produce the 
average phase space contraction — A„ if the segments beginning in ^(xq, n) produce A„. In 
analogy to Eq. (|52l) their steady state probability is given by: 

n-2 

7r,j(GA/"(xo),n) = /^iGM"(io) JJ^ P(^GAf"-'^(xo); ^ " ^ '^GM"^'^'^{xq)'-, k + 1) 

k=0 

X p{iGM{xo)] ^G(xJ, n) . (57) 

Again, this set of trajectories may be grouped together with the set of trajectories whose 
last step falls in the other region reachable from Gi^^, say i^,^ 7^ ir^^. The probability of the 
union a;(GM"(xQ), — 1) of these two sets takes the value 

n-2 

7rtj(GM"(xo),n-l) = /^iGA/"(xo) nP(^G'Af"-'={xo); ^ ~> ^GAf"-fc-i(i;o)' ^ + 1) 

fc=0 

— h^gmi^j^^x^^^^Pba Pda Pac Pcc Pbb Pdb Pad Pcd 

= l^GMi^j^.^^^^plAPccPBBPAD , (58) 



where we have made use, from Eq.(l36l). of the crucial relation icM'^ix^) = GMiM^-i^xg)^ with 
k = 1, ...,n, cf. Fig. 171 This contribution to the probability of producing — A„ in n steps 
mirrors the contribution to the probability of producing A„, given by Eq.(|56l). Taking the ratio 
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k=0 



k=l 



k=n-l k=n 




«Gil/"-i(ai) *GA/»(3i„) 



Figure 7. Sequence of visited regions in the forward (upper sequence) and time-reversed 
(lower sequence) dynamics. The lower sequence is determined by applying the composite 
map G o M onto the upper sequence, cf. Eq. ( |36] | (thick blue arrows). 



of these two contributions and writing the phase space contraction in terms of g units of size 
(j), cf. Eq.dSTl). one obtains 

7'a;(xo,n-l) _ /^J^^o VdA PbB PcC P AD 



with 



TT II rn^A ^nc ^no 

7ra;(GM"(xo),n-l) AiGMij^,,„_i(^^) PdaPccPbbPaD 

(59) 



am.n<a^= Pda^^Pbb^^Pcc^^Pad'''' < c^ma. , (60) 

where the upper and lower bounds ctmm and amax are {g, n) -independent positive numbers, 
which are found to be 



a. 



mm 4/ Ctmax ' (61) 

as we numerically tested by considering all possible values of corresponding to a trajectory 
segment visiting the regions i^^, 2Af"-i(xo), «M"(xo), for any i^^, im^^-^x^) and iM^ix^)- At 
the same time, Eqs. (l44l) and (|45l) imply the equality being at the heart of the A-FR, i.e. 

9 



Pbb 
Pec 

These results hold for all sets of trajectory segments starting in uj{xq, n — 1), related to their 
corresponding reversals starting in cj(GM"(xq), n — 1). Therefore, Eq.(|59]l holds as well 
for the total probabilities of producing A„ and — A„, because the ratio of the sums of the 
probabilities of the groups of trajectory segments producing A„ and — A„ equals the ratio of 
the probabilities of a single group, with corrections always bounded by and amax- 

To match this result with the A-FR Eq. (fT4)) . it now suffices to introduce the normalized 
quantity e„ = g(p/n{A) and to take the logarithm of the ratio of probabilities. 



en{x) — — < ——In— — — — — < en{x)+ — 7-r^-(62) 

n{A) n{A) fJ'{{x : en{x) e [—p — o, —p + 0)\ n{A) 
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In the n — oo limit, in which the allowed values of e„ become dense in the domain of the 
A-FR, one recovers the fluctuation theorem with p* = 4>/ (A) . 

6. Conclusions 

In this paper we have presented analytically tractable examples of dynamical systems in order 
to clarify some aspects of the applicability of the standard steady state fluctuation relation. In 
our case, there is no distinction between the so-called A-FR and il-FR, because the appropriate 
measure is the Lebesgue measure, in our case, cf. Eq.® E] dH. Our results show that the 
A-FR holds under less stringent conditions than those required by the Gallavotti-Cohen FT, 
which include time reversibility and existence of an SRB measure, i.e. a measure which 
is smooth along the unstable directions. This is of interest for applications, because strong 
requirements such as the Anosov property are hardly met by dynamics of physical interest, in 
general. 

To obtain this result, we have considered an example in which the involution representing 
the time reversal operator is discontinuous [fTTI and in which also the invariant measure is 
discontinuous along the unstable direction. Our discontinuities are mild, as discussed in the 
introduction, however, they illustrate how the validity of the A-FR may be extended beyond 
the standard constraints. Our proof capitalizes on the fact that the directions of stable and 
unstable manifolds are fixed and that the vertical variable does not affect the value of the phase 
space contraction rate. This fact has rather profound implications concerning the validity of 
the A-FR for cases in which time reversibility is more substantially violated. In fact, only 
the knowledge of the forward and reversed sequences of visited regions is required in order 
to verify the A-FR, rather than the more detailed knowledge of the forward and reversed 
trajectories in phase space. Thus, for instance, one easily realizes that our calculations may be 
carried out for a map of the form K = MoN, where M may refer to one of the maps Eqs.(fT5T) 
or (l33T l. while N does not contract or expand volumes and affects in some irreversible fashion 
the y-coordinate only. N can be constructed in several ways: For example, let M be the map 
Eq.(l33l). and assume that acts only on a vertical strip of width e in the region B, as follows: 



cf. FiglH] for a graphical representation. The map N is not reversible, according to the 
definition Eq.©; in fact, is not even a homeomorphism, as its inverse A^~^ is not defined, 
so neither is the inverse of the composite map K~^. Nevertheless, the A-FR still holds in 
this case, due to the existence of a milder notion of reversibility expressed by the relations 
Eq.(l36l). The latter entail that only a coarse-grained involution, mapping regions onto regions, 
is needed for the proof of the A-FR, rather than a local involution, mapping points into points 
in phase space, as defined by Eq.©. 




for X e [x,x + e] and y G 1]. 
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Figure 8. The map N defined in Eq.(l63]l, which spoils the reversibility of the model 
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